Linear and nonlinear waveguides induced by optical vortex solitons 
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We study, numerically and analytically, linear and nonlinear waveguides induced by optical vortex 
solitons in a Kerr medium. Both fundamental and first-order guided modes are analyzed, as well as 
the cases of effectively defocusing and focusing nonlinearity. 
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Optical vortex solitons are stationary nonlinear beams 
self-trapped in two spatial dimensions and carrying a 
phase singularity - a nonlinear analog of wave-front dis- 
locations . One of the possible applications of vortex 
solitons, as well as other types of spatial optical solitons, 
is to induce stable non-diffractive steerable waveguides 
that can guide and steer another beam, thus creating a 
reconfigurable all-optical circuit. 

In spite of many theoretical (see, e.g., and exper- 

imental H results on the soliton-induced waveguides in 
planar structures, waveguiding properties of optical soli- 
tons in a bulk medium have not yet been studied in detail. 
It is important to note that dark solitons are more at- 
tractive for waveguiding applications due to their greater 
stability and steerability jjj . Their two-dimensional gen- 
eralization, optical vortices, may have a number of poten- 
tial applications in trapping small particles [5] as well as 
waveguiding. Preliminary numerical results [3,0 indicate 
that some of the vortex waveguiding properties should be 
similar to those of planar dark solitons, however new fea- 
tures are expected. 

In this Letter we carry out a systematic analysis of the 
waveguides created by vortex solitons in a Kerr medium. 
We examine all the regimes. In the linear regime, the 
guided beam is weak and the analysis of waveguiding 
properties is possible by approximate analytical methods 
(see also Ref. R]). The nonlinear regime corresponds to 
the larger intensities of the guided beam and it gives rise 
to composite (or vector) solitons with a vortex compo- 
nent, identified and analyzed numerically. The effect of 
a focusing nonlinearity on a guided mode is also investi- 
gated. 

We consider the stationary evolution of two incoher- 
ently coupled beams with the frequencies u>i and lu 2 prop- 
agating in a bulk Kerr medium. Evolution of the slowly 
varying beam envelopes E\ and E 2 is described by the 
normalized equations ffl 
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of the beam with the frequency u>i j2 . The field ampli- 
tudes and the propagation distance z are measured in 
the units of /ei(no|"2|) 1/ ' 2 and feino, respectively, tiq is the 
linear refractive index, and n 2 characterizes an intensity- 
dependent change of the refractive index, An = n 2 |-E'| 2 - 

Equations ([!]) can describe two different physical sit- 
uations: the interaction of beams of two different fre- 
quencies (k and r\ are arbitrary), and the interaction of 
two orthogonally, linearly, or circularly polarized beams 
of the same frequency (77 = n = 1). In this Letter we 
consider the latter case. We assume that the beam E\ 
exhibits defocusing nonlinearity, whereas the beam E 2 
exhibits either defocusing (7 = +1) or focusing (7 = —1) 
nonlinearity, the latter may occur for two orthogonally 
polarized beams in a photorefractive medium |1. 

The parameter a > measures the relative strength 
of cross- and self-phase modulation effects. In gen- 
eral, its value depends on the physical nature of non- 
linearity, anisotropy of the medium, and polarization 
state of the beams. For example, for linearly po- 
larized beams in a Kerr-type material with a non- 
resonant electronic nonlinearity and cubic anisotropy, 
a = 2(2xxxyy + Xxyyx)/?<Xxxxx, Xijkl being the elements 
of the third-order susceptibility tensor. In isotropic ma- 
terials Xxxxx = 2>Xxxyy + Xxyyx and the values of a are 
fixed, namely a = 2/3 and a = 2 for linearly and cir- 
cularly polarized beams, respectively. Other values of a 
are possible, for instance in the case of nonlinearity due 
to molecular orientation a = 7 Jl0| . 

We assume that the envelope E\ carries a single- 
charged vortex on a nonvanishing background field \E\ \ = 
Eq, whereas the field E 2 is spatially localized, and it de- 
scribes either the fundamental or first-order mode guided 
by the vortex. Then, we look for radially symmetric so- 
lutions in the form: 



E 1 (R,^z)=r 1 1/2 E u{R)e~ lEo V" v , 
E 2 {R, <p; z) = V - 1/2 E a v(R)e- lWK)E ° 2z e mv , 



(2) 



where R = \J x 2 + y 2 and n, m are topological charges. 
Coupled equations for the real u and v become 
2 



A r u J u + u — [u + (a/ij)v ]u = 0, 



(3) 



where Aj_ is the transverse Laplacian, r\ = uj 2 /^ 2 , n = 
k 2 /k\, and k\^ 2 — 27r/Ai,2 is the propagation constant 
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-v + \v — j[v 2 + rjcru 2 ]v = 0, 
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where A r = d 2 /dr 2 + (l/r)d/dr is the radial part of 
the Laplacian, and r = RE$. For any a and A, Eqs. 
(||) admit a solution for the vortex soliton in the form: 
u = uo(r), v = 0. Here we are mainly interested in the 
case of a single charged vortex, thus n = ±1. For a fixed 
a and 77, the family of composite solitons (it, v) formed 
by a vortex and its guided mode is described by a single 
parameter A. For |u|/|w| = e <C 1, the second equation 
of (^) becomes a linear eigenvalue problem with an ef- 
fective potential created by the vortex. For 7 = +1, the 
potential is attractive, and it may support spatially lo- 
calized solutions - guided modes, each of them appearing 
at a certain cutoff value of A. For 7 = — 1, the potential 
is repulsive ( "anti-waveguiding case" ) , and no guiding is 
possible in the linear limit. 
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FIG. 1. Existence domain (shaded) for the fundamental 
mode guided by a vortex. Points a, b, c, and d correspond 
to the modal structure examples (a-d) shown below. Param- 
eters are: (a) - (c) a = 3, and (d) a = 0.8. Dash-dotted - the 
cut-off line calculated by variational approach. 

Away from the cut-off curve, the amplitude of the 
guided mode grows, and the linear waveguide approxi- 
mation is no longer valid. In this nonlinear waveguid- 
ing regime ||], the guided mode deforms the supporting 
vortex waveguide, and together they form a composite, 
vector soliton. To identify the families of such solitons, 
we find, numerically, localized solutions of (^) with the 
asymptotics: u — > 1, v — > as r — +00. As has been 
previously established for a = 2 there exist differ- 
ent solutions of this type, classified by the order of the 
guided mode. Here we describe the entire existence do- 
main for such solutions in the physically relevant region 
of parameters. 

Several examples of the coupled-mode solutions for 
7 = +1 and m = 0, i.e. when the vortex guides a fun- 
damental mode, are shown in Figs, l(a-d). The regions 
of existence of all possible solutions of this type are also 
presented in Fig. 1 on the parameter plane (tr, A) , where 
the points (a) to (d) correspond to the profiles below. 

In the parameter plane (A,er), the cutoff of the fun- 



damental mode is characterised by a curve <t(A) (thick 
curve, Fig. 1). Near this curve, the amplitude of the 
guided mode is small and the linear waveguiding regime 
is valid. Linearizing Eqs. (^|) around the vortex solution, 
we obtain two decoupled equations: 



A r uo 



7 u - u + u = 0, 



(4) 



A r v -v - Xv + W (r)v = 0, 



where only the terms ~ e are retained in the equation 
for the guided mode, and the effective guiding potential 
Wo(r) = — cruo(r) is defined by the vortex shape. Since 
the vortex profile uq (r) is not known in an explicit form, 
exact analytical solution of the eigenvalue problem for the 
mode v(r) is not available. However, an estimate for the 
guided-mode cutoff can be obtained by a standard vari- 
ational technique. First, we determine the vortex shape 
using a simple trial function uq(t) = tanh" (ar) and an 
effective Lagrangian for the first equation of ^) : 



duo 
dr 



-X 1 



(5) 



The averaged Lagrangian is then found as the integral 
(L) = a -1 Jq R L u (x)dx where x = ar, calculated in the 
limit R — > 00. Variation of (L u ) with respect to a yields 
the result: a = n 

_1 {/o°° x ^ ~ tanh 2 "(a;)] 2 dx}, valid for 
any n. This variational result provides the best fit to 
the numerical solution for n < 3. Then, the shape of the 
effective potential created by a single-charge vortex is de- 
termined by Uq(t) = tanh 2 (0.554r) wl - exp(— r 2 /4). 

In order to find the cutoff for the fundamental mode 
guided by the vortex from the second equation of (Q), 
we take into account the e 2 -order correction to the 
vortex-induced waveguide due to a nonzero v-amplitude 
W(r) — Wo + v 2 (r). Assuming a trial function of the 
form v(r) = B exp(— r 2 /b 2 ), and varying the averaged 
effective Lagrangian for the v-mode, 



+ Ar 



rv 
~2 



(6) 



with respect to b and B, we obtain a simple analytical 
expression for the cut-off: 
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1 + 24A + 16A 2 - (1 - 4A) X /1-4A) 2 



(7) 



A reasonably good agreement of this result with the nu- 
merically calculated cut-off can be seen in Fig. 1. 

The properties of the vortex-induced waveguides and 
the corresponding composite solitons are different de- 
pending on whether a < 2 or <j > 2 [see Fig. 1]. In the 
latter case, the amplitude of the guided mode grows with 
decreasing A [see Figs. l(a,b)] and the vortex waveguide 
becomes broader, as in the case of one-dimensional dark 
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solitons H . The amplitude of the guided mode reaches 
the vortex background far from the cutoff (dashed curve, 
Fig. 1) and goes above the background value before the 
solution expands and disappears (dotted curve, Fig. 1). 
Close to this limit the solution corresponds to a type of 
polarization domain walls [El . Its expansion with chang- 
ing A is very rapid, and hence the border of the existence 
region (dotted curve, Fig. 1) is determined by the accu- 
racy of numerical integration. For a < 2 the scenario is 
similar as A increases but the amplitude never reaches the 
background before the solution disappears on the other 
border (dashed curve, Fig. 1). 
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FIG. 2. Existence domain (shaded) and characteristic ex- 
amples for (a) a first-order mode guided by a vortex and (b) 
corresponding vector soliton. 
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FIG. 3. Existence domain of vector solitons (shaded) 

and two examples (a,b) of the vector solitons in the 

anti-waveguiding regime. 

Next, we examine the guiding of the first-order mode 
by an optical vortex, i.e. we take m = ±1 in Eqs. (|^). 
The existence domain for such solutions is shown in Fig. 
2 (left). For a = 1, Eqs. (|J) represent a two-dimensional 
generalization of a completely integrable one-dimensional 
Manakov system In this case, waveguiding prop- 

erties of a vortex resemble those of a one-dimensional 
Manakov dark soliton. Namely, the induced waveguide 
is single-mode, and the first-order mode exists only at 
A = a = 1. 

Although for a > 1 the vortex-induced waveguide is 
expected to guide higher-order modes, the corresponding 
composite solitons are likely to be dynamically unstable 
at least in the nonlinear regime shown in Fig. 2(b). As 
A decreases away from the cutoff, the amplitude of the 
guided mode grows and the waveguide becomes broader 
[see Fig. 2(b)] until its shape resembles that of a radi- 
ally symmetric domain wall B . The solution disappears 



when its localization region is no longer finite. 

At last, we briefly discuss the case of a focusing non- 
linearity, i.e. 7 = — 1 in Eqs. ([!]). Although no bound 
modes of the vortex anti- waveguide can exist in the linear 
regime (\v\ <C \u\), the numerical results suggest that for 
large intensities of the bright component mutual trapping 
is still possible. This effect is similar to nonlinear anti- 
waveguides assisted by an external potential p2| . Sta- 
tionary solutions shown in Fig. 3 (right) for a fixed a, 
appear when an absolute value of A (A < for 7 = — 1) 
exceeds a certain threshold value (solid curve, Fig. 3). 
Away from this border, the amplitude of the bright soli- 
ton increases. 

Our results can be readily generalized to the case of 
beams of different colours, r\ ^ 1. For the limit of lin- 
ear waveguiding, the effective coupling constant changes 
to 770- [see Eqs. (^)], and the properties of the vortex 
waveguide become dependent on the frequency ratio. 

In conclusion, we have analysed the waveguiding prop- 
erties of optical vortex solitons in Kerr media and shown 
that they crucially depend on the relative strength of the 
cross- and self-phase modulation effects. We have identi- 
fied families of composite solitons formed by a vortex and 
its guided mode with or without a topological charge and 
introduced a simple technique of determining the cutoff 
of the guided mode by a variational approach which can 
be applied to other problems of nonlinear waveguiding. 
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